Electroproduction of the S , n(1535) Resonance at High Momentum Transfer 
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The differential cross section for the process p(e, e'p)n has 
been measured at Q 2 — 2.4 and 3.6 (GeV/c) 2 at center-of- 
mass energies encompassing the 5*11(1535) resonance. The 
latter point is the highest-Q 2 exclusive measurement of this 
process to date. The resonance width and the helicity-1/2 
transition amplitude are extracted from the data, and evi- 
dence for the possible onset of scaling in this reaction is shown. 
A lower bound of » 0.45 is placed on the 5n(1535) — > pn 
branching fraction. 

PACS Numbers: 13.60.Rj, 13.40.Gp, 13.60.Le, 14.20.Gk 



tainty in the Sii(1535) transition amplitude is limited 
by knowledge of the full width and branching fraction to 
the 77. We report here on a measurement of the reaction 
e + p — ► e' + 5n(1535) — * e! + p + r\ and an extraction 
of the helicity-conserving transition amplitude A^, 2 at 

Q 2 = 2.4 and 3.6 (GeV/c) 2 . We also use a recent anal- 
ysis of inclusive (e, e') data to put a lower bound on the 
<Sii(1535) — > pr\ branching fraction. 



II. THE EXPERIMENT 
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I. INTRODUCTION 

Baryon electroproduction allows the measurement of 
transition form factors, which test models of hadronic 
structure in ways that static baryon properties alone 
cannot. Recently much effort has gone into attempts 
to reproduce observed transition form factors over a 
large range of four-momentum transfer. At low four- 
momentum transfer, or Q 2 , the focus has been on in- 
corporating relativistic effects into the constituent quark 
model (CQM) Q, using light-front g-f§ and other ffl 
approaches. At higher Q 2 , perturbative QCD (pQCD) 
sum rule calculations j?J and valence pQCD || have been 
employed. The applicable range in Q 2 for these various 
approaches is not clear. 

Among the most interesting of baryon case studies 
is the <Sii(1535) resonance, which is one of the most 
strongly excited states over all Q 2 , and which is eas- 
ily isolated because it is the only resonance that has a 
large branching fraction to the rj. The reproduction of 
the Sii(1535) form factor has become a goal of many 
models, but the effort has been hampered by a lack of 
precise electroproduction data. In addition, the uncer- 



The experiment was performed in Hall C of the 
Thomas Jefferson National Accelerator Facility (Jeffer- 
son Lab), shown in Figure [|. The Short Orbit Spectrom- 
eter (SOS), which is a resistive QDD device, was used 
to detect electrons. The High Momentum Spectrometer 
(HMS), which is a superconducting QQQD spectrom- 
eter, was used to detect protons. Figure || shows the 
HMS detectors, which include drift chambers (DC1 and 
DC2) for determining track information, scintillator ar- 
rays (S1X/Y and S2X/Y) for triggering and time-of-flight 
measurement, and a threshold gas Cerenkov and elec- 
tromagnetic calorimeter for particle identification (P1D). 
The SOS detectors are configured similarly. 

The incident electrons had energies E = 3.245 and 
4.045 GeV for the Q 2 = 2.4 and 3.6 (GeV/c) 2 points, 
respectively. At each of the two Q 2 points, the electron 
spectrometer was fixed in angle and momentum, thus 
defining a central three-momentum transfer q and direc- 
tion of a boosted decay cone of protons. The proton 
spectrometer was stepped in angle and in momentum to 
capture as much of this decay cone as possible. Data 
were obtained at 33 (21) kinematic settings at the low 
(high) Q 2 point. 

Target protons were provided in the form of liquid 
hydrogen at 19 K flowing through a target of length 



4.36 cm. The relative current of the electron beam was 
measured by two resonant-cavity current monitors, which 
were calibrated periodically using the absolute beam cur- 
rent measurement of a parametric current transformer. 
The combined measurement had an absolute accuracy of 
(7 = 1.5%. 

Electrons were identified in the SOS using the 
Cerenkov detector and lead-glass calorimeter (see Fig- 
ure ||). In the HMS, protons were separated from pions 
using the time of flight measured between two pairs of 
scintillator arrays (see Figure^). In both spectrometers, 
tracking information was obtained from the drift cham- 
bers. Details of the experiment and analysis are given in 
Ref. Q , and information on a simultaneous measurement 
of the A(1232) can be found in Ref. @. 

III. DATA ANALYSIS 

The data were corrected for trigger and PID inefficien- 
cies (< 1%), track reconstruction inefficiencies (« 5%), 
computer and electronic dead times (< 5%), current- 
dependent target density changes (~ 3%), and protons 
undetected due to interactions in the detector stack 
(~ 3%). The data were binned in W, cos 8* , (j>* , and 
M 2 (with 6, 10, 6, and 20 bins, respectively). Here W 
is the invariant mass; 9* is the polar angle between the 
direction of the rj and the three-momentum transfer q 
in the center-of-mass (cm.) of the resonance; <fi* is the 
azimuthal angle of the rj with respect to the electron scat- 
tering plane; and M 2 is the square of the missing mass 
for p(e, e'p)X. The rj mesons were identified in the fi- 
nal state using M 2 . Figure JH] shows the missing mass 
distribution for a typical kinematic setting. 

Modest backgrounds in M 2 due to accidentals (« 2%, 
shown in Figure |]) and protons penetrating the HMS 
collimator and magnet apertures (« 4%) were measured 
and subtracted from the data. The remaining contin- 
uum background in missing mass was due to multi-pion 
(nir) production (ranging from 30% to 50% of the reso- 
nance data) and a small (< 2%) contribution from target- 
window interactions. Two independent techniques were 
used to subtract these remaining background events. The 
first technique fitted a polynomial plus peak in M 2 to the 
data in each (cos#*, <f)*) bin (integrated over the W ac- 
ceptance for that kinematic setting), and then subtracted 
the background contribution from each bin. The sec- 
ond technique scaled a Monte Carlo-generated nir back- 
ground to match the data above and below the missing- 
mass peak and then subtracted this background from 
each (W, cos#*, </)*) bin. 

Three different models were used to simulate the nir 
background in the Monte Carlo: ep — > e' pir + n~ , ep — > 
e' A ++ 7r~ — * e' pn + 7r~ , and a crude approximation of 
three-body phase space. The Monte Carlo simulation was 
also used to simulate multiple scattering and ionization 
energy loss, and to correct for experimental acceptance 



and the effect of radiative processes. Once the nir back- 
ground was subtracted from both experimental and sim- 
ulated spectra, the experimental yields were corrected to 
account for finite Q 2 acceptance. The differential cross 
section was then given by the ratio of experimental to 
simulated yield in each (W, cos#*, </>*) bin, normalized 
by the simulation resonance cross section for that bin. 

The cross sections obtained using the different nix mod- 
els and the two background subtraction techniques all 
agreed within 2%; both the following figures and our fi- 
nal results were obtained using the first subtraction tech- 
nique together with a background generated by combin- 
ing two of the nn models. Figure |^ shows data and fits for 
several typical (cos#* , (j)*) bins of one kinematic setting. 
Figure shows the result of fits for several kinematic set- 
tings, where for each setting we have integrated both the 
data and their respective fits over the sixty individual 
(cos6>* , 4>*) bins. 

Using similar techniques we verified the well-known 
1 H(e, e'p) cross section |nj to within 2%. 

IV. RESULTS 

The five-fold differential cross section for the ep — > 
e'pn process may be expressed as the product of the 
transverse virtual photon flux T T (Hand convention [ fl2|| ) 
and the cm. cross section for the electroproduction of 
the prj pair: 

da da 

Previous data indicate that the cm. ^f v p — > pn 
cross section is dominated by S-waves arising from the 
5n(1535) ]T^jr^| . This dominance was confirmed by the 
present data, which showed that terms other than S'-wave 
were less than 7% and consistent with zero within the sta- 
tistical uncertainty of the data. Angular distributions for 
the Q 2 = 3.6 (GeV/c) 2 data are shown in Figure ||. 

From S'-wave fits to the angular distributions, the to- 
tal cross section was calculated (at each Q 2 point) as a 
function of W: 

otot (W) = 4tt ( lv p -> pn) . (2) 

This cross section, which consists of resonant and nonres- 
onant parts, was fitted with a relativistic Breit-Wigner 
plus nonresonant background curve, 

Otat(W0 = cr rcs(H / ) + a m (W) 

2 \ P ;\w w 2 r 2 R 

rcs m p K (W 2 -W^) 2 + W 2 i T 2 (W) 

+ B m - W thl , (3) 

where Wr is the resonance mass, is the full width, 
A 2 es and B nT are the Q 2 -dependent magnitudes of the 



resonant and nonresonant terms, K is the equivalent real 
photon energy [K = (W 2 — m 2 )/(2m p )], and p* is the 
three-momentum of the r\ in the cm. of the prj system. 
The pi] production threshold is at Wthi ~ 1486 MeV 
(in the lowest W bin). At both values of Q 2 , the 
fitted value of the phenomenological nonresonant term 
(B nT y/W — Wthr) was consistent with zero (with an un- 
certainty of 1% of the resonant term). 

The energy-dependent resonance width T(W) of Eq. || 
was parameterized in terms of the branching frac- 
tions b n (= r, ; /r fl at Wr), b^, and according to 



Walker |15|. At present the Particle Data Group (PDG) 
gives an estimated range for the 77 branching fraction of 
0.30 < b n < 0.55 @. Therefore, fits to a ICS {W) were 
made assuming three sets of values for the branching 
fractions (b^ : b^ : b^), which we define as Fits 1-3, 
respectively: (0.55 : 0.35 : 0.10), (0.45 : 0.45 : 0.10), and 
(0.35 : 0.55 : 0.10). A consequence of the pr\ threshold 
is that the fit to o~ rcs (W) cannot constrain the branching 
fractions Q (i.e., the three fits result in curves that are 
virtually indistinguishable). 

Based on a branching fraction constraint presented be- 
low, we consider Fit 1 (b v = 0.55) to a 1 - cs (W) to be the 
preferred fit. The fits for both Q 2 points are shown in 
Figure |[ With the Fit 1 branching fractions, we obtain 
a full width T R = (154 ± 20) MeV. This width changed 
less than 10 MeV over the range of branching fraction 
assumptions. The uncertainty is statistical added in 
quadrature with systematic. Our result agrees with the 
PDG estimate (« 150 MeV) |16|, and appears lower than 
the recent Mainz measurement, = (203±35) MeV 
(see inset of Figure |fj|). These recent results disagree 
with the value of Tr — (68 ± 7) MeV obtained from 
the high-Q 2 measurement of Ref. pj. The form of the 
Breit-Wigner parameterization used by the three groups 
is essentially the same, and so does not account for the 
differences in Tr. 

As noted above, the fit to <r ICS (W) cannot constrain 
the branching fraction b n , but a comparison between this 
work and a recent fit to inclusive (e, e') scattering |l7| ] 
can. The fit by Keppel et al. models the inclusive cross 
section in terms of transverse resonant (0T res ) and non- 
resonant {<JT nr ) contributions using 



da 



dflrdEi 



T T [a TnT (1 + eRnr) + <J Tr , 



(4) 



In that work, the resonant contribution from each of the 
three resonance regions (assumed to be entirely trans- 
verse) is fit using a Breit-Wigner form. The transverse 
component of the nonresonant contribution is fit using 
the phenomenological form 



v = Y J C n (Q 2 )(W- W thl ) n -^ 



(5) 



where the C n {Q 2 ) are fourth-order polynomials in Q 2 . 
The longitudinal component of the nonresonant cross sec- 



tion, which enters through the longitudinal-to-transverse 
ratio R nr , is taken from a fit to deep inelastic data flg|| . 

The resonant part of the second resonance region is 
dominated at low Q 2 by the Di3(1520). At higher 
Q 2 , however, the Sii(1535) begins to dominate, and by 
Q 2 = 4 (GeV/c) 2 it is expected that the (1535) is 
responsible for over 90% of the resonant cross section at 
W w 1535 MeV juj. Assuming that the resonant part 
of the inclusive cross section is the incoherent sum of the 
resonant contributions of the various decay channels, we 
can use the inclusive and exclusive resonant cross sections 
to put a lower bound on b^ ||: 



b v > 



s (5ji ->pr]) 
,8 (inclusive) 



(6) 



where both cross sections are taken at W ss 1535 MeV. A 
value of b n = 0.55 results in good agreement between the 
high-Q 2 point of this work and the inclusive fit; a value 
of b n = 0.35, on the other hand, implies an inclusive cross 
section 50% greater than the fit to the measured inclusive 
cross section, which is strong evidence that the branching 
fraction is not this low. With the incoherent summation 
ansatz given above, and assigning a 10% uncertainty to 
the inclusive fit, we find a lower bound of b n = 0.45 with a 
95% confidence level. Assuming complete S±i dominance 
at Q 2 = 4.0 (GeV/c) 2 , we find a best fit of b v = 0.55. 

Neglecting resonances other than the Sii(1535), we re- 
late the amplitude A^ 2 to cr ros by [ ]l3| , |l9| 



Al /2 {Q 2 ) = 



w r t r a res (g 2 , W R ) 



2 m p brj 



1 + eR 



1/2 



(7) 



Here e is the longitudinal polarization of the virtual pho- 
ton, and R — a L /a T . For R we assumed a parame- 
terization based on a quark- model calculation p(J. The 
expected impact of the longitudinal-to-transverse ratio R 
on the final physics result is small: a 100% error in the 
assumed value [~ 4% at Q 2 — 2.4 (GeV/c) 2 ] corresponds 
to an uncertainty of less than 1% in the quoted value of 

A l/2 ■ 

Table |I| gives final results for Fits 1-3. The uncertain- 
ties are systematic and statistical added in quadrature; 
for Ay 2 we included estimates for the uncertainties from 
Tr and b v , which were obtained by varying these quanti- 
ties over reasonable ranges (150-200 MeV and 0.45-0.6, 
respectively) and studying the effect on the helicity am- 
plitude. 

Table |l| lists the dominant sources of systematic un- 
certainty in the measurement and their impact on the 
differential cross section and on the helicity amplitude. 
The uncertainty in i s given as a range, where the 

largest uncertainties are for the highest W bins. 

Figure [l(] shows the helicity amplitude results, along 
with points calculated from previous e p — > e'p rj data and 
some theoretical predictions. All data points in the figure 
were calculated using Eq. assuming Tr — 154 MeV 



TABLE I. Results. The uncertainties are systematic (in- 
cluding estimated uncertainties in Tr and b v for A®, 2 ) and 
statistical added in quadrature. The top A® , 2 result is for 
Q 2 = 2.4 (GeV/c) 2 , the bottom is for Q 2 = 3.6 (GeV/c) 2 . 
Fit 1 is preferred for reasons discussed in the text. The 'best 
value' for b v assumes Su dominance at Q 2 = 4 (GeV/c) 2 . 



Quantity 



Fit 1 



Fit 2 



Fit 3 



W R [MeV] 
T R [MeV] 



^1/2 



[io- 
[io- 



GeV 
GeV 

r,/r« 



-1/2] 
-1/2, 



1532 ± 5 
154 ± 20 

50 ±7 

35 ±5 
> 0.45; 



1527 ± 5 
150 ± 19 
55 ±8 
39 ±6 
best value 



1521 ± 5 
147 ± 19 
63 ±9 
44 ±6 
0.55 



TABLE II. Dominant sources of systematic uncertainty, 
not including Tr and b v (which affect A^ 2 ). 



Quantity 


Fractional uncertainty (a) in 

d(l* ^1/2 


Monte Carlo nn model 


1% to 7% 


1% 


nir subtraction 


1% to 6% 


1% 


Knowledge of E 


1% to 10% 


0.8% 


Knowledge of 8 e 


0.2% to 11% 


1% 


Experimental Acceptance 


1% to 6% 


1% 



and b v = 0.55; if either assumption is wrong, all data 
points will scale together. Not included for any of the 
data points in the figure are the uncertainties in T# and 
b v . Note the good agreement between the high-Q 2 point 
of the present work and the inclusive fit for b v = 0.55; 
assumption of a lower branching fraction shifts the data 
up relative to the inclusive fit. 

The present result differs from previous work in both 
the strength and the slope of the 5n(1535) form factor; 
most notably, we find a cross section 30% lower than that 
of Ref . |uj [found by interpolating the results of this work 
to Q 2 = 3 (GeV/c) 2 ]. This difference is reduced in the 
amplitude by the square root relating A^ 2 to the cross 
section (Eq. |^). Although the present data were taken at 
a different value of e than those of Ref. p4|, a longitudinal 
cross section is not responsible for the difference between 
the two measurements; a value of R « 2.3 (which is ruled 
out at low Q 2 [^1 22 1) would be necessary to account for 
the discrepancy. 

Of the various CQM curves shown in Figure \tu, none 



exhibit a slope as shallow as that of the data. Those that 
indicate an amplitude at Q 2 ~ 3 (GeV/c) 2 roughly con- 
sistent with experimental data also predict excess ampli- 
tude at lower Q 2 . Our data also have consequences for a 
recent coupled-channel model for the Sn(1535) Hj; the 
proposed quasi-bound KH (five quark) state is expected 
to have a form factor that decreases more rapidly than is 
observed. 

Figure |l| shows the quantity Q 3 Ay 2 for the 5n(1535), 
which is predicted by pQCD to asymptotically approach 
a constant at high Q 2 M. As has been pointed out else- 



where [ p8| , such scaling might be due to non-perturbative 
contributions. While there is no strong scaling evident 
in the figure, our data indicate that Q 3 A^ 2 may be ap- 
proaching a constant value by Q 2 ~ 5 (GeV/c) 2 . 



V. CONCLUSIONS 

We have presented the results of a precise, high statis- 
tics measurement of the ep — > e'pr) process at W ~ 
1535 MeV and at Q 2 = 2.4 and 3.6 (GeV/c) 2 |§. The 
contribution of terms other than ,5- wave multipoles is ob- 
served to be less than 7%, which is consistent with previ- 
ous measurements. More importantly, the cross section 
obtained from the new data is about 30% lower and in- 
dicates a full width twice that of the only other exclusive 
measurement at comparable Q 2 [[uj. 

While the new data exhibit no strong perturbative sig- 
nature, they do have a Q 2 dependence that is markedly 
different than the older high-Q 2 measurement. Even 
given the new (lower) cross section obtained from this 
measurement, however, relativized versions of the quark 
model fail to reproduce the Q 2 dependence seen experi- 
mentally. 

A comparison of the new high-Q 2 datum (the highest 
in existence) with a recent inclusive analysis indicates an 
Sii(1535) — > prj branching fraction of at least b v = 0.45 . 
Using b v — 0.55 we obtain Tr = 154 MeV and a new 
measurement of A^ 2 (Q 2 ) (see Table |). 
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FIG. 1. A plan view of the Hall C end station at Jefferson 
Lab. The electron beam enters from the left, and the scatter- 
ing takes place in the cryogenic target placed in the beamline. 
In this experiment, outgoing particles were detected by two 
magnetic spectrometers: the Short-Orbit Spectrometer (SOS) 
was used to detect electrons and the High-Momentum Spec- 
trometer (HMS) was used to detect protons. 

FIG. 2. A side view of the HMS detector stack, as seen 
from the door of the detector hut. The detected particles 
travel from left to right (along positive z). 

FIG. 3. The response of the SOS calorimeter and Ceren- 
kov for events of a typical data run. The calorimeter response 
_B ca i is the total energy deposited normalized to the particle 
momentum, while the Cerenkov response N p . c . is the number 
of photo-electrons detected. The events at N p . c . — are n~ 
(note the peak at E CZj \ « 0.25). The events at N p . c . > 0, 
-Ecai > 0.7 are electrons. The events at N p . c . > 0, -E ca i ~ 0.3 
are caused by tt~ that produced knock-on electrons that trig- 
gered the Cerenkov. Note that the z axis is on a log scale. 

FIG. 4. Velocity from time of flight (/3hms) and coincidence 
time (the difference in time of arrival for the two spectrome- 
ters) for events of a typical data run. The band of events at 
Phms ~ 1 are 7v + , while those at /?hms ~ 0.8 are protons. The 
real proton coincidences are at t = ns, and the nominal 2 ns 
radio frequency structure of the beam is visible in the adja- 
cent accidental peaks. The 1ow-/3hms tail emanating from the 
real coincidence peak is most likely due to protons undergoing 
interactions in the detectors after the drift chambers. 

FIG. 5. A plot of Ml for one kinematic setting. The peak 
at M 2 « 0.3 (GeV/c 2 ) 2 corresponds to undetected r\ mesons 
in the final state (the peak at M 2 « 0.02 (GeV/c 2 ) 2 corre- 
sponds to 7T°, the subject of Ref. ^)). Note the presence 
of the multi-pion background as well as the radiative tail ex- 
tending to the right of the rj peak. 

FIG. 6. Fits to the M 2 distribution for several typical 
(cos#* , </>*) bins, one kinematic setting. 

FIG. 7. Results of background fits for several typical kine- 
matic settings. Data are on the left and the corresponding 
Monte Carlo result is on the right. Each figure shows the 
integration of sixty individual (cos#^ , <f>^) bins and their re- 
spective fits (like those shown in Figure ^|). The solid line 
is the sum of the background and peak fits; the dashed line 
shows the background only. The lines at the bottom of the 
data plots show the small contribution from the accidental 
coincidence and HMS collimator backgrounds. 

FIG. 8. Angular distributions for the Q 2 = 3.6 (GeV/c) 2 
data. Each plot shows the cos 6^ distribution for a single 
(W, (j>^) bin. The rows correspond to different bins in W , 
the columns to different bins in <f)^ . Data corresponding to 
<(>*! = ±90 degrees are not shown; the out-of-plane experimen- 
tal coverage was complete only for the lowest W bin (where 
the data looked similar to that in the (j>^ bins shown here), and 
was almost nonexistent at higher W . The lines are S-w&ve 
fits to the data. 



FIG. 9. Fit 1 to <7 r cs(W) for the two Q 2 points of this work 
(errors on the data are statistical only). Note the presence 
of the pn threshold. The inset shows the VF-dependence of 
this cross section as measured by the present work (solid line, 
F R = 154 MeV), Ref. (l3) (dashed line, T R = 203 MeV), and 
Ref. H (dotted line, V R = 68 MeV). The curves in the inset 
have been normalized to the same magnitude. 



FIG. 10. The helicity amplitude A\ /2 {Q 2 ) of the Sii(1535), 
measured via ep — » e'pr/, together with some theoretical pre- 
dictions. The data points ( ll3l|l4|,E]l-G5[ and the present 



work) were calculated using T R = 154 MeV, b v = 0.55, and 
the parameterization of R referenced in the text. The errors 
shown on previous data are statistical only. The errors shown 
for the present work include both statistical and systematic 
uncertainties, with the exceptions noted in the text. The the- 
oretical curves of Refs. [0-pU2q] are based on variants of the 
CQM. The curve from Ref. is the result of a pQCD cal- 
culation. The curve from Ref. |l7| is a fit to inclusive data. 



FIG. 11. The quantity Q S A\ /2 {Q 2 ) for the Si^^S). 
The dot-dashed line is an exponential fit to the cross 
section given by the two points of the present work 
(<r rcs = 16.5 exp [-0.565 Q 2 ] fih, where Q 2 is in [(GeV/cfj), 
and the solid line is a fit to inclusive data (as in Figure fLOh . 
The errors that are shown, and the assumed values for T R , 
b„ , and 7?, are the same as in Figure 
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